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' Recent experiments on Mn-doped ZnS nanocrystals have shown unusual magnetization properties. 

We describe a nearest-neighbor Heisenberg exchange model for calculating the magnetization ratios 
^ of these antiferromagnetically doped crystals, in which the dopant atoms are distributed inhomoge- 

neously within the nanocrystal. This simple inhomogeneous doping model is capable of reproducing 
the experimental results, and suggests that interior dopant atoms are localized within the crystal. 
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, Transition metal doping is now possible for a variety of nanocrystalline semiconductor materials, e.g., CdSf^*^ 
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I. INTRODUCTION 



CdSe)^, ZnSe^ and ZnS&i In bulk dilute magnetic semiconductors (DMS), Heisenberg coupling of magnetic dopants, 
G such as Mn, occurs via anion-mediated superexchange^^ but it is not completely understood how this is modulated 
^ ' by finite nanocrystal (NC) size effectsi^ 

o : 

^ ' In a recent experiment, Tsuji et al. synthesized Mn-doped ZnS NCs and determined the ratio between the magneti- 



zation at a finite field and the expected magnetization for the given number of dopants per NC, defined as M^T/Mgat^ 



■ The initially prepared NCs (designated "as prepared") had a 1.2% Mn dopant concentration, an average of 4.3 Mn 
dopants per NC (mostly on the surface), and M^T/Mgat = 91%. Washing the sample with acid removed the surface 
atoms, yielding a sample (designated "HCl-washed" ) with a 0.10% Mn concentration, an average of 0.3-0.5 Mn dopants 
per NC (presumably on the interior), and M^T/Msat — 77%. The effect of antiferromagnetic coupling of the Mn ions 
is to reduce the magnetization ratio. However, the "as prepared" NCs, with their higher dopant concentration (and 
hence increased chance of exchange interaction between dopants), show a higher M^r/Mgat than the "HCl-washed" 
samples, with their order of magnitude lower dopant concentration. Tsuji et al. posit that Mn ions on the surface 
do not interact magnetically. However, in this paper we show that even with uniform Heisenberg exchange constants 
for surface and interior atoms, taking into account the nonhomogeneous distribution of dopant atoms within the NC 
alone is sufficient to explain the magnetization results. 
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II. MODEL 

The experiments of Tsuji et al. are conducted with 3 nm diameter ZnS NCs^i Using the zinc-blende model of Lippens 
and Lanoo, this corresponds to a 729 atom NC, with alternating "shells" of anion and cation atoms, summarized in 
Table Qpfi Assuming random incorporation of dopant atoms during the growth process, the distribution of dopant 
atoms within the crystals obeys a binomial law as noted by Counio et al.^ where the probability, p{Nd), that a NC 
contains Nd dopant atoms, is given by 

piNd)^{ (1) 



Nd 



where 



\ni/ 

is a binomial coefficient, x is the mean dopant concentration, and ngu^.s is the number of cation sites into which the 
dopant atom may substitute itself. This is easily generalized to consider separate "inner" and "outer" volumes within 
the NC, consisting of riinner and Uouter sites and having dopant concentrations Xmner and Xouter, respectively. This 
is illustrated schematically in Figure ^ 

The magnetization ratio, M^x/Msat (using the notation of Tsuji et al), is related to the observed susceptibility, 
Xobs and that expected of N independent dopants of susceptibility xii and can be broken down into contributions of 
various dopant cluster geometries by the expression 

XobsfJ-oB 



{M^T/Msat) 



NdXif^oB 



where i runs over the possible cluster types. Pi is the probability of occurence of a given cluster of dopant atoms, 
and Xi — Xi/xi is the normalized susceptibility of cluster-type i, as calculated by Kreitman et al}^- Xi is a function 
of temperature [T = 5K in the experiment) and the nearest-neighbor Heisenberg exchange constant, which we 
take to be — J/fcfc ~ lOK for Mn in bulk ZnS.^'^ With (anion-mediatcd) nearest-neighbor interactions, one obtains 
from a statistical mechanical treatment of the spin configurations as a function of —k}jT/J the normalized cluster 
suceptibilities for all clusters up to three atoms as: Xsingie = 1-0 (by definition), Xpair ~ 0.0, X open triangle = 1-0, 
and X closed triangle — 0.1."'^^ A morc detailed description is given in the Appendix. The values of the Xi are essentially 
constant for 5K < — J/kjy < 50K, so our results are independent of the exact value of J. 
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For the bulk case, the probabihties of occurence for the various clusters have been derived by Behringerii^ For 
nanocrystalline systems, this is unsatisfactory, due to the finite size and shape effects on the probability of forming 
various clusters; recently, Suyver et al. have performed numerical simulations and derived analytic results for the 
probability of dopant pairs in spherical NCs assuming homogeneous doping. Less than 1% of Mn atoms are part 
of three-dopant clusters for the concentrations in our calculations, so the pair-only truncation would make a minor 
contribution to the magnetization ratio. However, we have included them in our calculations, as the determination 
is not particularly difficult. More important is the possibility of an inhomogeneous distribution of dopant atoms 
within the NC, illustrated schematically in Figure Q We performed numerical simulations, in which Eq.(^ was 
used to determine a number of dopant atoms to be randomly placed within a given set of cation shells. P,; was 
evaluated by counting all the types of (anion-mediated) nearest-neighbor clusters in a statistical sampling of 10^ 
dopant configurations. Then Eq.(|3Jl was used to determine the magnetization ratio for the ensemble of doped NCs. 

III. RESULTS 

Using the model of a 729 atom, 3.0 nm ZnS NC, we varied the local dopant concentration, keeping the concentration 
over the total NC equivalent to the magnetically determined concentration, Xcurie, from experiment. The shells are 
explicitly described in Tabled and the results are shown in Table HTl The higher average number of dopant atoms, Nd, 
in our calculations reflects the fact that our model NC is slightly larger than that assumed by Tsuji et al; otherwise 
we would expect x curie to underestimate the number of dopants by a few percent^ 

We initially performed calculations in which the dopant atoms were uniformly distributed over all 420 cation sites; 
this results in a higher magnetization ratio for the low-concentration samples than for the high-concentration samples, 
contradicting the experimental results. However, we understand the acid washing procedure as removing any dopant 
atoms on the outermost shell of the NC. The x^c reported for the "HCl- washed" experimental samples corresponds 
to the number of dopant atoms in the interior divided by the total number of cation sites, ntotaU in the NC (including 
the surface cation sites). But in the "HCl- washed" case, the remaining dopants are not in the outermost shell, but 
are constrained to be in the riinner', we can determine Xinner using the relation xjsicntotai = Xinnernmner- Similar 
stoichiometric reasoning allows us to calculate Xouter ■ The calculations in Table ^ show the results of constraining 
the inner sites to consist of the one, two, three, and four cation shells, both with and without doping into the exterior 
(fifth) cation shell. For cases with an outer shell and higher xnc (modelling the "as prepared" experimental samples), 
Xouter > Xinner rcduccs the probability of pairing between surface and interior dopant atoms, increasing M^x/Mgat as 
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compared to the uniformly distributed case. However, the magnetization ratio for these high xj^ic samples is relatively 
insensitive to any restriction of the sites available for the interior dopant atoms, since the number of dopant atoms 
in this exterior shell is much higher than the number of dopant atoms in the interior, and to a first approximation is 
similar to the result obtaining by having all the dopant atoms on the exterior (since Xouternouter ^ Xmnernmner)- 

Constraining the interior dopant atoms to reside within some locally concentrated sphere (for the purpose of our 
model, concentric with the NC), is qualitatively sufficient to achieve the proper effect of reducing the magnetization 
ratio for the washed samples with respect to the unwashed samples (as seen by comparing the respective rows of the 
"With Exterior Shell" and "Without Exterior Shell" calculations in Table ITl)l . In particular, constraining the interior 
dopant atoms to reside on the innermost cation shell nearly reproduces the experimental magnetization ratio results 
for both the "as prepared" (93.4% theory, 91% experiment) and the "HCl-washed" (80.1% theory, 77% experiment), 
despite the simplicity of the model. 

IV. CONCLUSION 

The inhomogeneous dopant distribution model reproduces the magnetization behavior observed by Tsuji et 
aZ.without assuming different magnetic interactions for surface and interior dopant atoms. As noted in Section 
im at the experimental temperature of 5K, Heisenberg exchange constants in the range of 5K < — J/fcf, < 50K are 
indistinguishableii^ Further study on the effects of dopant geometry and NC size on J may rely on electronic structure 
calculations, such as the recently developed ZILSH method of O'Brien and Davidsonii& 

Our results suggests that dopant atoms in the interior of the NC are localized. It has been noted previously that 
MnxCdi_xSe NCs appear to effectively zone-refine out the dopant atoms, since the energetic barrier to rearranging 
the lattice is not as insurmountable as in a bulk crystal.'^ This suggest that there may exist a region near the center of 
the NC in which the energetic barrier is high enough to prevent migration of dopant atoms to the surface, surrounded 
by a region in which dopants migrate to the surface during synthesis and processing, as shown schematically in Figure 
^ This may be further investigated by molecular dynamics simulations. 
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VI. APPENDIX 



Below is a brief summary of the cluster susceptibility model, following Kreitman et alr^ and Smarfeii. The Heisen- 
berg Hamiltonian for nearest-neighbor interactions is 



H=-2J S,-S,+g^iBHoY,S- 

nn{i,j) i 



(4) 



where J is the exchange constant, nn{i,j) indicates that only nearest-neighbor i and j pairs are included, Sj is the 
spin- vector for site i, g is the g-factor, /i^ is the Bohr magneton, Hq is the applied field, and is the z component of 
the spin on site i. There are three types of nearest-neighbor clusters containing three atoms or less: pairs (-D), closed 
triangles {E) and open triangles (iV). The eigenenergies of Eq.Q) for these configuratations are respectively, 

Ed = -J[JiJ + l)-2SiS+l)]+g^iBHoM 

< J < 25, |Af I < J 
Ee = ~J[ji3 + l)-2SiS+l)]+g^iBHoM 

0<k<2S,\k-S\<j <k + S, \M\ < j 
En = -J[jij + l)-k{k+l)-SiS + l)]+giiBHoM 

0<k<2S,\k-S\<j <k + S, \M\ < j, 

using the standard notation of spin angular momentum. For n € {D, E. N^, the free energy, F, is given by. 



(5) 



I3F ^ In^e" 



(6) 



where (3 = I/ZcbT, /cb is the Boltzmann constant and T is the temperature. The susceptibility, x, is given by 

X 



dHh H„=o 



Ho=0 



Pigf^B? E„ (E, Sife-P^- ( g^^B E„ (E. s^fe-^^- 



E„ e 



Hn=0 



J Ho=0 



(7) 
(8) 

(9) 

(10) 
(11) 
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The susceptibility of a set of noninteracting (J = 0) spins, xs, is simply, 

xs^lpigf^sfsis + i). (12) 

We can then use this to definte a normalized susceptibility X = xlxs, ^^nd likewise normalized susceptibilities of each 
cluster type i, as Xi = Xi/xsi used in Eq.©- 
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FIG. 1: Schematics of the inhomogeneous dopant distribution model proposed in this paper. Xi„rier and Xouter are the dopant 
concentrations for the inner and outer sheUs respectively, riinner and Uouter are the number of cation sites enclosed within 
the inner and outer shells, respectively. Black circles show a random configuration of dopant atoms. For the low dopant 
concentration ("HCl washed") samples, the restriction of doping to a small inner shell increases the probability of pairs as 
compared to the case where the dopant atoms are uniformly distributed over the entire NC, hence reducing M^t / Mgat- For 
high dopant concentration ("as prepared") samples, Xouter > Xinner reduces the probability of pairing between surface and 
interior dopant atoms as compared with the uniformly distributed case, hence increasing Mz,t / Msat ■ 
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TABLE I: Description of the zinc-blende NC model of Lippens and Lanoo— for a 3.0-nm diameter ZnS crystal. 
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Uniformly Distributed 

420 1.20 0.00 1.20 5.03 90.3 

420 0.10 0.00 0.10 0.42 99.1 
With Exterior Shell 

216 0.19 204 2.27 1.20 5.02 93.7 

92 0.46 204 2.27 1.19 5.01 94.5 

28 1.50 204 2.27 1.19 5.00 94.0 

4 10.48 204 2.27 1.17 4.91 93.4 
Without Exterior Shell 

216 0.19 204 0.00 0.10 0.42 98.2 

92 0.46 204 0.00 0.10 0.41 96.2 

28 1.50 204 0.00 0.10 0.40 90.6 

4 10.48 204 0.00 0.07 0.31 80.1 
Exterior Shell Only 

0.00 204 2.27 1.09 4.60 94.3 
All Dopants on Exterior 

0.00 204 2.47 1.20 5.02 93.8 

0.00 204 0.21 0.10 0.42 99.5 
Experimental Results- 

As prepared 1.2 4.3 91 

HCl-washed 0.10 0.3-0.5 77 

TABLE II: Calculated results, for lO'' NC sample. Definitions of the columns are given in Sections 2 and 3. We found the 
values of M^t / Maat and Nd to vary only slighlty as compared to calculations using 10'* NCs. 



